Elaboration of the a-Model Derived from the BCS Theory of Superconductivity 
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The single-band a-niodel of superconductivity [H. Padamsee, et al., J. Low Temp. Phys. 12, 387 
(1973)] is a popular model that was adapted from the single-band Bardeen-Cooper-Schrieffer (BCS) 
theory of superconductivity mainly to allow fits to electronic heat capacity versus temperature T 
data that deviate from the BCS prediction. The model assumes that the normalized superconducting 
order parameter A(T)/A(0) and therefore the normalized London penetration depth Al(T')/Al(0) 
are the same as in BCS theory, calculated using the BCS value obcs ~ 1.764 of a = A(0)/fcBT'c, 
where fee is Boltzmann's constant and Tc is the superconducting transition temperature. On the 
other hand, to calculate the electronic free energy, entropy, heat capacity and thermodynamic critical 
field versus T, the a-model takes a to be an adjustable parameter. Here we write the BCS equations 
and limiting behaviors for the superconducting state thermodynamic properties explicitly in terms 
of Q, as needed for calculations within the a-model, and present plots of the results versus T and 
a that are compared with the respective BCS predictions. Mechanisms such as gap anisotropy 
and strong coupling that can cause deviations of the thermodynamics from the BCS predictions, 
especially the heat capacity jump at Tc, are considered. Extensions of the Q-model that have 
appeared in the literature such as the two-band model are also discussed. Tables of values of 
A(T')/A(0), the normalized London parameter A(T)/A(0) and AL(r)/AL(0) calculated from the 
BCS theory using a — qbcs are provided, which are the same in the Q-model by assumption. 
Tables of values of the entropy, heat capacity and thermodynamic critical field versus T for seven 
values of a including OfBCS are also presented. 

PACS numbers: 74.20. De, 74.25. Bt, 74.20. Rp, 74.20.Fg 



I. INTRODUCTION 

The 1957 Bardeen-Cooper-Schrieffer (BCS) micro- 
scopic single-band clean-limit mean-field theory of super- 
conductivity describes superconductivity as arising from 
an indirect attractive interaction between Cooper pairs of 
electrons mediated by the electron-phonon interaction.^ 
In the BCS theory, a Cooper pair has spin 5* = (a 
spin singlet) with zero angular momentum correspond- 
ing to s-wave superconductivity. The BCS theory is a 
weak-couphng theory, which means that 






< 1 and 



A(0) 
fceOo 



«1, 



(1) 



where Tc is the superconducting transition temperature, 
0D is the Debye temperature, fce^D is the maximum 
phonon energy within the Debye theory,^ /cb is Boltz- 
mann's constant, and A(r) is the superconducting order 
parameter versus temperature T, which is the activation 
energy for single quasiparticle (electron and hole) exci- 
tations of the superconducting ground state and is often 
just called the superconducting gap. The actual super- 
conducting energy gap 2 A is centered on the Fermi en- 
ergy of the metal and is thus qualitatively different from 
the energy gap in a semiconductor where the energy gap 
is between the top of the valence band and bottom of the 
conduction band. The BCS theory makes precise pre- 
dictions of A(T) and the magnetic penetration depths 
derived from it, together with predictions of thermo- 
dynamic quantities which include the superconducting 
state electronic free energy Fes(T), thermodynamic criti- 
cal field Hc{T), entropy Scs{T) and heat capacity Ccs{T). 



The BCS theory predicts that A(T) decreases contin- 
uously and monotonically to zero as T^ is approached 
from below, resulting in a second-order phase transi- 
tion at Tc- A finite discontinuous increase ("jump") 
ACe(Te) = C{Tc) - Cc„(Te) in Cc(r) occurs on enter- 
ing the superconducting state from the normal state with 
decreasing T. The jump is given in the weak-coupling 
limit by ACc/'jnTc ~ 1-43 for all BCS superconductors 
because it is a law of corresponding states, where the 
normal-state electronic heat capacity is Cen{T) = jnT 
and 7n is called the Sommerfeld electronic heat capac- 
ity coefficient. This value of AC'e/jnTc was confirmed 
for some superconductors such as Al and Ga, but other 
superconductors showed ACo/7nT'c > 1.43 such as the 
value 2.7 for Pb.'^ The enhanced jumps were subsequently 
determined to arise in superconductors that arc not in 
the weak-coupling limit, termed moderate- or strong- 
coupling superconductors.^ There exists no rigorous the- 
oretical expression for CcsiT) for such superconductors 
that can be routinely used to fit experimental heat ca- 
pacity data such as the heat capacity jump at Tc. 

To provide a model for fitting experimental Ccs(T) 
data such as ACo(Tc)/7nT'c and other thermodynamic 
properties for moderate- and strong-coupling supercon- 
ductors, Padamsee, Neighbor and Shiffman introduced 
the single-band a-modcl in 1973,^ which is able to 
fit CcsiT) data not only for strong-coupling supercon- 
ductors with ACc(Tc)/7nT'c > 1-43 which was the 
original motivation, but also for superconductors with 
ACc(Tc)/7nTc < 1.43 which can arise from anisotropy of 
A in momentum space (see, e.g., Ref. 6, and Sec. VIII 
below). The a-model assumes that the normalized gap 



A(T)/A(0) and therefore the normahzed London pene- 
tration depth AL(r)/AL(0) are the same as in the BCS 
theory which are calculated using the BCS value aecs ~ 
1.764 of a = A(0)/A;Brc- On the other hand, the nor- 
malized superconducting state Scs{T)/jnTc, Ccs{T)/jTc, 
Fcs{T)/jT^ and Hc{T)/.J^ are calculated taking a 
to be an adjustable parameter in the corresponding BCS 
equations to allow fits to experimental data that deviate 
from the BCS predictions. 

Thus the a- model is not self-consistent, but provides a 
popular model with which experimentalists can fit their 
electronic superconducting state thermodynamic data 
that deviate from the BCS predictions and to quantify 
those deviations. However, few of the many papers re- 
porting use of the a-model to fit experimental Cca{T) 
data explain how the theoretical values for the fits were 
obtained. Indeed, the original exposition by Padem- 
see et al. in Rcf. 5 is not completely developed. The 
main purpose of the present paper is to write the BCS 
equations for the thermodynamic properties explicitly in 
terms of a so that these equations can be directly used 
to numerically calculate and plot these properties versus 
both T and a. We also discuss predictions of the BCS 
theory itself for comparison with the predictions of the 
a-model. Tables of values of the various superconducting 
state properties for both the BCS theory versus T and 
the a-model versus T and a are provided. 

In order to introduce and explain calculations needed 
for the a-modcl, the predictions of the BCS theory^'^ 
and their limiting behaviors at low and high tempera- 
tures are first discussed. These BCS superconducting 
state predictions are described in Sees. II-VI including 
plots of calculated data. The equations are written in 
terms of the a parameter of the a-model, which are then 
used in Sec. VII to calculate ^es, Ces, and He versus T 
and a, which are compared with the BCS theory pre- 
dictions. The role of gap anisotropy as a mechanism for 
producing reduced heat capacity jumps ACe/7n21; < 1-43 
is discussed and calculated for several example cases in 
Sec. VIII. Extensions of the single-band a-model that 
have appeared in the literature are briefly discussed in 
Sec. IX. A summary is given in Sec. X. In this paper, 
we use Gaussian cgs units throughout.^ 

In the Appendix, tables are provided of values 
of A(i)/A(0), the London parameter A(i)/A(0) and 
AL(i)/AL(0) calculated from the BCS theory using a = 
Q^BCS I where t = T jT^^ which are the same in the a-model 
by assumption, and a table of the BCS prediction of the 
Pippard penetration depth Ap(f)/Ap(0) versus t. Also 
provided in the Appendix are tables of the normalized 
values Scs(t)/jnTc, Cos(t)/7nT'c and iJc(i)/V7n^ ver- 
sus t for seven values of a including aecs ■ These tables 
supplement the very popular table in Ref. 8 of supercon- 
ducting state properties versus temperature predicted by 
the BCS theory. In the areas of overlap, the values in our 
tables agree with those in Ref. 8.^ 



II. BCS GAP EQUATION AND Tc 



The BCS gap equation is 



de , f E 
—- tanh 



E \2kBTj N{0)V' 



where 



E ^ Ve^ + A2 



(2a) 



(2b) 



is the energy of an electron excited above the supercon- 
ducting energy gap and of a simultaneously excited hole 
below the gap with respect to the Fermi energy Ep at 
the center of the gap, both of which are defined to be 
positive, e is the normal-state single-particle energy, V 
is the elcctron-phonon coupling constant and A is as- 
sumed to be the same everywhere on the Fermi surface 
(isotropic s-wave superconductivity). The excited elec- 
trons and holes are together termed quasiparticles. It 
is assumed that kBOu/EY <C 1, where Ep is the Fermi 
energy measured from the bottom of the electron con- 
duction band. In the BCS theory, the energies e and E 
arc measured with respect to Ep, so Ep = for those 
energies. A^(0) = N{Ep) is the normal-state electronic 
density of states at e = Ep = for a single spin direc- 
tion, which is assumed to be constant over the energy 
range of interest in the theory. The full energy gap that 
develops below the superconducting transition tempera- 
ture Tc is 2A(T), with Ep in the middle of the gap. Thus 
A is the activation energy for a quasiparticlc excitation 
and 2A is the energy required to break up a Cooper pair 
via simultaneous excitations of an electron and a hole 
quasiparticlc. The energy E in Eq. (2b) approaches the 
normal state energy e for e 3> A. 

To determine Tc, one sets A = and T = Tcin Eqs. (2) 
and takes the weak-coupling limit Tc/6d ^ 1, yielding 
the BCS equation for Tc given by 



Tc = ( 1 Od cxp - 



1.134 Ooexp - 



iV(0)F 
1 



(3) 



iv(o)Fy ' 



where 7e ~ 0.5772 is Eulcr's constant. Defining the sym- 
bol 



"BCS = -^ ~ 1.764, 



Eq. (3) can be written 



1 , ^ 290 

m 



A^(0)1/ VaBCsTc 



(4) 



(5) 



Inserting this into the right side of Eq. (2a) gives a second 
form of the gap equation 



feee 



— tanh I — — — I = In I — ■ 

E \2kBT) \aBcsTc 



(6) 



III. BCS SUPERCONDUCTING ORDER 
PARAMETER A(T) 

Taking the hmit T — >■ and the weak-coupHng hmit 
Tc/Q-D — >■ in Eq. (6) allows A(0) to be evaluated as 



A(0) 
ksTc 



aBcs- 



(7a) 



The full superconducting energy gap at T = is therefore 



2A(0) 
kv.Tr, 



2aBCS 



27r 
e7E 



3.528. 



(7b) 



Inserting the expression for olbcs in Eq. (7a) into the 
right side of the gap equation (6) gives a third form of 
the gap equation 



L 



''"'^^ de ( E 

-— tanh — — — 
E \2kBT 



In 



2fcBeD 



A(0) 



(8) 



In order to calculate the temperature dependence of A 
we define dimensionless reduced variables 



A 



A 
A(0)' 



A(0)' 



T 



From Eq. (2b), the excited quasiparticles have reduced 
energy 



E 



E 

A(0) 



vi2 + A2. 



(10) 



Using Eq. (7a) and the definitions (9) and (10), the gap 
equation (8) expressed in dimensionless variables is 



^^^''""""itanhf^Uln 



2fcBeD 



A(0) 



(11) 



where now the ratio /i;B0D/A(O) occurs both in the upper 
limit to the integral on the left side and in the argument 
of the logarithm on the right side. 

Numerically solving Eqs. (10) and (11) for A at fixed 
values of t in the asymptotic weak-coupling regime with 
fcBOD/A(0) = 500 > 1, where the calculated A(i) is nu- 
merically independent of the precise value of fcB0D/A(O), 
yields the A versus t data shown in Fig. 1. A list of A{t) 
values is given in Table II in the Appendix, where a log- 
arithmic scale in 1 — t is used to provide a high density of 
A values for t -^ 1 where A varies rapidly. Also shown 
in Fig. 1 is a plot of A^(i) versus t, from which one sees 
that the order parameter near Tc is A(t ^ 1) ~ ^/l — t, 
a t dependence characteristic of mean-field behavior as 
further discussed below. 

The low-t behavior of A{t) can be obtained from 
Eqs. (10) and (11) by doing an integration by parts and 
then utilizing the weak-coupling limit A:b0d/A(O) — )• oo. 
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FIG. 1: (Color online) Normalized BCS energy gap 
A(r)/A(0) and its square versus reduced temperature T/Tc 
as obtained by numerically solving Eq. (11) in the weak- 
coupling limit (1). 



yielding 
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(12) 



from which expressions for the temperature derivatives 
d"A{t — s- 0)/(it" can also be calculated. 

An expression for computing dA^/dt for < t < 1, 
which we need below to calculate Ccs{t), is obtained by 
taking the t derivative of the gap equation (11) and solv- 
ing for dA^/dt, yielding 



dA^t) 



J^ sech^ {g)de 



dt 

where 


Jo 


tsccli^ig) t2ta,nli(g) 


de 




aBCS\/e^ + A2 
n . 





(13a) 



2t 



(13b) 



Alternatively, one can use Eqs. (40) and (42b) below to 
calculate dA^{t)/dt, which give the same numerical re- 
sults as Eqs. (13). 

For T -> Tc, one obtains 



dA^it) 



dt 



t->i 



'7C(3) 



-3.016 



(14) 



where ^(x) is the Riemann zeta function. From Fig. 1 
and Eq. (14), for < — >■ 1 one has 



A^{t^l)^A{l-t). 



(15) 



Equations (14) and (15) give 

A{t^ 1) = x/A^/T"^ 



(16a) 



1.737 VT~~^, 



which is a t dependence of the order parameter char- 
acteristic of mean-field theories as noted above. Equa- 
tion (16a) can also be written as 



A(^->1) 
knTr 



aBcsA{t 



7C(3) 



Vi^ 



3.063 ^/]~~t, 



(16b) 



where we used the expression for aecs in Eq. (4). For 
later use, we combine Eqs. (14) and (15) to obtain 



t dA^ 



A2 dt 



1 



1-t 



(17) 



Using Eqs. (13), dA^/dt and dA/dt at t ^ are cal- 
culated to be 



(18a) 



dA2(f^0) 

dt 

dA{t -^ 0) 

dt 


2V27raBCS aecs/t (, 

1 dA^it^O) IdA^it^O) 
^ 2A dt ^2 dt 




\/27raBCS aecs/t fi 



t3/2 



(18b) 



These two results are the same as obtained directly from 
Eq. (12) by taking the respective derivative and utiliz- 
ing only the leading term ^ \/t in the prefactor of the 
exponential. 
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FIG. 2: (Color online) Normalized electronic entropy 
Sc{T)/^nTc versus reduced temperature t = T/Tc in both the 
superconducting (T/Tc < 1) and normal {T/Tc > 1) states 
for three values of a = A(0)/fcBT'c including the BCS value 
a = Qbcs in Eq. (4). The superconducting state data are cal- 
culated using Eq. (20a). The normal-state behavior at T > Tc 
is given by Scn/'ynTc — T/Tc and the extrapolation to lower 
T is the long-dashed line. 



where the Fermi-Dirac distribution function is (with 
Sf = 0) 



f = fiaBcs,E,t) 
with 



E ctBCsE 



(20c) 



k^T 



t 



IV. BCS ELECTRONIC ENTROPY AND HEAT 
CAPACITY 

The normal-state Sommerfeld electronic specific heat 
coefficient 7„ is related to A^(0) by^ 



2^2fc2 



7n 



iV(0). 



(19) 



The 7n is the measured value and therefore both 7,1 
and iV(0) contain the same enhancements by many-body 
electron-electron correlations and the electron-phonon in- 
teraction. Using Eq. (19), the superconducting-state 
electronic entropy Sea and heat capacity Cos are ex- 
pressed in terms of 7112!; and reduced variables according 
to 



Sea{t) ^ 6q;|cs 
7nTc 



ttH 



fiaBcs,E,t)[E 



E 



dl (20a) 



Ccs(0 „ 6a 

7nT'c 



3 

BCS 



ttH 



;,!-/) |!-1^ I..-, 



t 2 dt 



(20b) 



The normalized electronic entropy S^/jnTc versus 
reduced temperature t calculated numerically using 
Eq. (20a) is plotted as the solid black curve in Fig. 2. 
One sees that the superconducting and normal state en- 
tropies are the same at Tc, which demonstrates that the 
superconducting transition is second order with no latent 
heat. The superconducting state entropy is lower than 
the normal state entropy for < t < 1, showing that the 
superconducting state is more ordered than the normal 
state in this temperature range. The superconducting- 
state entropy approaches zero exponentially for t — > 
due to the superconducting energy gap for quasiparticle 
excitations. 

To obtain Cos at a particular t one first determines 
A at that t using Eq. (11), inserts that into Eqs. (13) 
to determine dA'^/dt at that t, and then inserts these 
two quantities into Eq. (20b) and does the integral there. 
All integrals are done numerically. A plot of the BCS 
prediction of Co versus t for a = aecs is shown as the 
black solid curve in Fig. 3(a). The heat capacity jump 
ACo(To) on cooling below Tc is given by Eq. (20b) as 



ACo(ro) 

711^0 



3aBcs^\ dA 



27r2 



dt 



(21) 
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FIG. 3: (Color online) Calculations using the a-model. 
(a) Normalized electronic heat capacity Cc{T)/'y-nTc versus 
temperature T in both the superconducting (T/Tc < 1) and 
normal {T/T^ > 1) states for three values of a = A(0)/A;Brc 
including the BCS value a = aacs in Eq. (4). The super- 
conducting state data are calculated using Eq. (20b). The 
normal-state behavior at T > Tc is given by Ccn/fnTc = T/Tc 
and the extrapolation to lower T is the long-dashed line, (b) 
Dependence of the specific heat jump ACe(rt;)/7nTc at Tc 
on a. The BCS value aecs ~ 1.764 in Eq. (4) is indicated. 



Substituting the expressions for aecs from Eq. (4) and 
dA'^/dt\t=i from Eq. (14) into (21) gives 



Aa(re) _ 12 

7nTc 7C(3) 



1.426. 



(22) 



At low temperatures t < 0.2 where A « 1 (see Table II 
in the Appendix), the heat capacity is given by BCS as 

C.s{t) 3«lcs [3^^(^^^g/,) ^ K,{aBcs/t)], (23) 



Tn^'c 



27r2f2 



where K^ (x) is the modified Bessel function of the sec- 
ond kind. This heat capacity arises from excitations of 



electron and hole quasiparticles, and at these low temper- 
atures does not include a contribution from a change with 
temperature of the superconducting condensation energy 
^ A^. We have verified that numerical data generated 
using Eq. (20b) are in precise agreement with Eq. (23) 
at i < 0.2. The expansion of Eq. (23) at low t is 



Ces(i) 3V2aBcs /"bcsA^/^ 



"fnTc 



7r3/2 

11 

1 + ^ 



t 



225 / t 



"BCS/ 128 V"BCS 



(24) 



X e 



-aBCs/t 



Thus at low t, Ces decreases exponentially to zero with 
decreasing t. Using aBcs/t = A(0)/fcBr, this T depen- 
dence is seen to arise from excitations of electron and 
hole quasiparticles above and below the superconducting 
energy gap, respectively, with activation energy A(0) for 
both types of quasiparticle. The numerical prefactor is 
3A/2aBcs/7i-3/^ ~ 1.344.10 

BCS fitted their calculations of Ccs(i)/7nTc for f < 0.7 
using Ec^. (20b) by the expression Ces(t)/7nTc = ae"''/* 
and obtained a — 8.5 and b = 1.44. Because a was 
a constant independent of t and the fit was not done 
in the low-i limit, the fitted value b = 1.44 in the ex- 
ponent is not equal to the low-i limit value aBCS = 
A(0)/fcBT'c w 1.764. The fit was done to compare their 
theoretical prediction with experimental data that were 
not in the \ow-t limit. Within BCS theory, one could evi- 
dently obtain A(0) from fits of experimental Cos(T) data 
for T < 0.7 Tc by the expression Cos(T) = a'e"'''/''''^ 
using A(0) « 1.22 6', where 1.22 = 1.76/1.44. 



V. BCS ELECTRONIC FREE ENERGY F^s 
AND THERMODYNAMIC CRITICAL FIELD H^ 

The electronic Helmholtz free energy F^ is^^ 

Fe = C/e - TS„ (25) 

where Ue is the electronic internal energy. The normal- 
state free energy is 

+2 

(26) 



Fen{t) 



t' 



InT^ 2 • 

The free energy i^og of a BCS superconductor in the weak- 
coupling limit is 

Fos{t) 



InT^ 



(27) 
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For A — > at t — > 1 one regains the normal-state free 
energy, i.e., 

Te,(t=l) _ Fenft=l) _ 1 
7nT2 ^^^2 2- 



(28) 
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FIG. 4: (Color online) Thermodynamic critical field He cal- 
culated using Eqs. (26), (27) and (29), normalized by its zero- 
temperature value Hc{0) in Eq. (30), (a) versus reduced tem- 
perature t = T/Tc and (b) versus t^ , according to the a- model 
for the a values listed, including the BCS value a = aecs in 
Eq. (4). 



The superconducting state has a lower free energy than 
the normal state for < t < 1, demonstrating that the 
superconducting state is the ground state in this temper- 
ature interval. 

The thermodynamic critical field Hc{T) is defined 
in terms of the free energy difference at zero applied 
magnetic field between the normal and superconducting 
states versus temperature as 






fcn(t) - fcs(t) 
7nT2 • 



(29) 



Since H^ is expressed in cgs units of Oe, where 
1 Oe^ = 1 erg/cm'^, the free energy is normalized to unit 
volume with units of erg/cm'^ and 7n is expressed in units 
of erg/cm'^ K^. The value of Hc{T = 0) obtained from 



Eqs. (26), (27) and (29) is 



gc(0) /6 /— _^^ 

r-rr = \ - aecs = V 67r e ^^ w 2.438, 



(7nT<2) 



(30) 



where the expression for aecs in Eq. (4) was used to 
obtain the second equality. A plot of Hc/Hc{0) versus t 
calculated numerically using Eqs. (26), (27) and (29) for 
a = CKBCS is shown in Fig. 4(a) and versus t^ in Fig. 4(b). 
From the latter figure one sees that Hc{t)/Hc{0) « 1 — f^ 
as noted by BCS. A fist of Hc{t)/Hc{0) values versus t 
for a = ctBCS is given in Table V in the Appendix.^ 

One can also determine Hc{t) from the difference be- 
tween the normal- and superconducting-state entropies 
Scn{t) — Scs{t) according to^^ 



7nT<? 



Stt 



Stt 



t>cn{t ) Jcs{t ) 



InTc -fnTc 



dt' 



Scsjt') 



dt' 



(31) 



where Scn{t')/^nTc = t' , and Scs{t')/jnTc is given in 
Eq. (20a). Equation (31) is used to extract Hc[T) from 
experimental heat capacity data after subtracting the lat- 
tice and possible magnetic contributions. With in BCS 
theory, we find that the normalized iJc/v7n^c versus t 
calculated from Eq. (31) is identical to that calculated 
from the free energy in Eq. (29) and plotted in Fig. 4(a) 
for a = aBCS) as must necessarily be the case. 



VI. MAGNETIC FIELD PENETRATION 
DEPTH 

The magnetic field penetration depth \{T) of a super- 
conductor is defined as the length scale of penetration 
of an external magnetic induction into a semi-infinite su- 
perconductor with the field applied parallel to the flat 
surface of the superconductor. ^'^^ Taking the external 
held direction as the z-direction and the direction per- 
pendicular to the surface as the x-direction, the general 
definition is^'^ 



A(T) = 



1 



B, 



CXt 2 J 



B^{x,T)dx, 



(32) 



where Hext = Boxt are the applied magnetic field and 
magnetic induction and B is the magnetic induction in- 
side the superconductor. Two limiting regimes for A are 
discussed by BCS. The London limit with A = Al corre- 
sponds to Co/Al <C 1, where 



Co 



hvf 



1 



hvY 



7rA(0) Traecs kB,Tc 



0.1805- 



hv-p 

kv.Tr 



(33) 



is the BCS coherence length, which is the minimum 
length scale over which A can change significantly at low 
temperatures t <^1, and v-p is the Fermi velocity (speed). 



Here local electrodynamics is used, which apply to type- 
II superconductors. Most superconducting compounds 
are in this regime. The other limit is the Pippard limit 
^0/-^? 3> 1 with Ap = A in which nonlocal electrodynam- 
ics is important and which applies to extreme type-I clean 
superconductors such as pure Al with Tc = 1.2 K."^ In the 
BCS paper, quasiparticle scattering by impurities is not 
included in any of the calculations. This "clean limit" 
corresponds to -^ :^ max(^o, A), where i is the mean free 
path for quasiparticle scattering by impurities. BCS com- 
mented that A is expected to increase with decreasing £, 
as elaborated by Tinkham.^ 



A. London Local Electrodynamics 

In the London local limit of the electrodyamics where 
^0/''^ ^ Ij using the notation in Eq. (32) the London 
equations yield^ 



B,(x,T) = Bext.(a; = 0)e- 



-s/Al(T) 



(34) 



where Al = A is the London penetration depth. The 
superconducting current density J(r) is related to the 
magnetic vector potential A(r) in the London gauge by 



J(r) 



1 



rA(r), 



cA(t) 
where for a free electron gas one has 



A(0) = 



TO 



7T,e^ 



(35) 



(36) 



Al(0) = 



A(0)c2 
47r 



Amne^ 



-, (37) 



Wn 



inne'^ 



(38) 



c is the speed of light in vacuum, m is the electron mass, 
n is the normal-state conduction electron density, e is the 
fundamental electric charge and ujp is the plasma angu- 
lar frequency of the conduction electrons in the normal 
state. 2 

The BCS solution for the London parameter A(i) and 
the corresponding supcrfluid density ps (t) in the two- fluid 
model is 



Psjt) _ A(0) 
Ps(0) A(i) 



= 1 - m, 



where 



1 f°° 
I{t) ^ -J sech" {z)dy, 



with z 



y^ + {aBcs^/ty 



(39) 

(40a) 
(40b) 



In the limit i — >■ 0, one can set A = 1 and replace sech(2;) 
by 2e~^ in the integral (40a). Then the integral can be 
evaluated analytically, yielding 



I{t ^ 0) 



27ra 



t 



BCS -QBCs/t 



(41) 



BCS relate A(t) to the temperature derivative of the 
gap according to 



m 



1 



t dK{t) 



1 



t 



A(0) ^ A(i) dt ' 2A2(t) 
Then using Eq. (39) one obtains 



dK^jt) 
dt 



t dA{t) 



t dA2(i) 



= 1 - 



1 



(42a) 



(42b) 



A(t) dt 2A2(i) dt ' l-I{t) 

In the limit t ^- 1, combining Eqs. (17) and (42b) gives 
1 1 



= 1 



1- I{t^l} ' 2(1 -i) 

The BCS prediction for AL(i) is, using Eq. (39), ^^ 



(43) 



m 



1 



Al(0) V A(o) yT^T(t)^ 

which can be written 

AAL(t) _ Al(<)-Al(0) _ 1 



Al(0) 



Al(0) 



V^^W) 



-1. 



(44a) 



(44b) 



From Eqs. (42a) and (44a) and Fig. 1, \i,{t) / \i,{Q) is 
equal to unity at i = and diverges to cxj for i — > 1. 
BCS commented that superconductors that are not in 
the London limit at low temperatures would be expected 
to come into that limit for i —J- 1 as A diverges. 

At low temperatures t <C 1 one has I{t) ^ 1 
from Eq. (41), so one can Taylor expand Eq. (44a) 
about I{t) = and use Eq. (41) to obtain the low-i 
approximations ^ ^ 

AL(t-J'O) _, , 1^,^, _, , /TTCkBCS 



Al(0) 



= 1 + ^m = 1 



-aBCs/t 



2i 



(45a) 



and 



For t - 
[l-I{t 



AAL(t->0) ^ / ^aBCS ^.-OBcs/t 
Al(0) \j 2t ^ 



(45b) 



1, inserting the expression in Eq. (43) for 
1)]^-*^ into (44a) gives 



Al(< 



Al(0) 



2(1 -t) 



(45c) 



which diverges to 00 at i = 1 as noted above. Very 
close to Tc, Eq. (45c) becomes \i^{t -^ 1)/Al(0) w 
1/^/2(1^. 
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FIG. 5: (Color online) (a) Variation of the BCS London pene- 
tration depth ^ ^SJ = j^qy — 1 versus reduced temperature 
t — T/Tc in the London limit ^o/Al ^ 1 computed using 
Eqs. (40) and (44). The low-T approximation in Eq. (45a) 
and the high-T approximation in Eq. (45c) are also shown. 

+2 



(b) Inverse penetration depth squared, )^ rj. , versus t, 
i^ and f*, as indicated. 



A plot of AAL(t)/AL(0) versus t from Eq. (44b) and 
a nmnerical solution of Eqs. (40) is shown in Fig. 5(a), 
together with the low- and high-T limiting behaviors in 
Eqs. (45a) and (45c), respectively. A list of values of 
the normalized London parameter I{t) = 1 — -Mr from 
Eqs. (39) and (40) and normalized London penetration 



depth 



AAL(t) 
A(0) 



AL(t) 
Al(0) 



1 



- 1 from Eqs. (40) 



and (44a) versus t is given in Table VI in the Appendix. 

To test power law temperature dependences, in 
Fig. 5(b) are plotted [Al(0)/Al(<)]^ versus t, t"^, t^ and f^. 
None of the four plots are linear over the whole tempera- 
ture range, but the t'^ dependence comes closest overall.^ 



FIG. 6: (Color online) Plots as in Fig. 5, but derived from 
Eqs. (46a)-(48) for the Pippard penetration depth Ap in the 
Pippard limit ^o/Ap 3> 1. 



B. Pippard Nonlocal Electrodynamics 

In the Pippard limit ^o/X ^ 1 for which nonlocal elec- 
trodynamics is appropriate in clean extreme type-I su- 
perconductors, the BCS prediction for the t dependence 
of A = Ap for diffuse scattering of quasiparticles from the 
surface of the superconductor is 



Ap(0) 
Ap(i) 



A(i)tanh 



aBcsA(^) 
2t 



1/3 



where 



Ap(0) 
Al(0) 



2^Al(0) 



1/3 



(46a) 



(46b) 



Since ^q/X > 1 in the Pippard limit, Eq. (46b) 
gives Ap(0) > Al(0), as also expected from gen- 
eral arguments.^ A list of values of AAp(i)/Ap(0) = 
[Ap(t)/Ap(0)] — 1 is given in Table VII in the Appendix. 



For i — >■ 0, we use the lowest-order expansion of 
A in Eq. (12) and the large argument approximation 
tanh(a:) = 1 — 2e~'^^ in Eq. (46a) to obtain 



Ap(0) 
Ap(t->0) 



> 



2TTt 
CkBCS 



-OBCs/t 



(47) 



For t ^ 1 where A — ^ 0, one can set tanh(x) = x in 
Eq. (46a) and use Eqs. (4), (14) and (15) to obtain 



Ap(0) 



Ap(f 



1 



47re'^'= 

7c(3y 



(1-i) 



1/3 



1.3856 (1-t) 



1/3 



(48) 

A plot of AAp(i)/Ap(0) versus t calculated from 
Eq. (46a) is shown in Fig. 6(a) along with the low- and 
high-T approximations and their extrapolations using 
Eqs. (47) and (48), respectively. Plots of [Ap(0)/Ap(i)]^ 
versus i" {n = 1-4) are shown in Fig. 6(b). As noted by 
BCS, the calculations as displayed in Fig. 6(b) are ap- 
proximately described by the Gorter-Casimir two-fluid 
model with n = 4 given by 



m 

m 



i-t\ 



(49) 



A comparison of Figs. 5(a) and 6(a) shows that Ap(i — > 
1) diverges more slowly than Xj^it — )■ 1). However, as 
discussed above, a superconductor in the Pippard limit 
^o/A ^ 1 at low t would cross over to the London limit 
^o/A <C 1 as Tc is approached due to the divergence of A 
at T^. 



VII. PREDICTIONS OF THE «-MODEL 

The superconducting state entropy, heat capacity, free 
energy and thermodynamic critical field are calculated 
within the a-model by replacing aecs in Eqs. (20), (21), 
(23), (24), (27) and (30) by an adjustable parameter 
a = A(0)/fcBT'c. As noted in the introduction, the BCS 
t dependence of the reduced gap A(t) in Fig. 1 calcu- 
lated using Eq. (11) is retained, which is carried out 
using the BCS value aecs ~ 1.764 in Eq. (4). From 
Eqs. (42a) and (44a), the normalized London penetra- 
tion depth AL(i)/AL(0) is uniquely related to A(i) and is 
hence also the same in the a-model as in the BCS theory. 

The electronic entropy So normalized by ^nTc deter- 
mined using Eq. (20a) is plotted versus t in Fig. 2 for 
three values of a including the value a = aecs • One sees 
that the superconducting- and normal-state entropies are 
the same at Tc for a = 1.5 and 2, as is the case for the 
BCS value aecs ~ 1.764, so the transition at T^ remains 
second-order within the a-model. The electronic specific 
heat CchnTc obtained from Eq. (20b) is plotted versus 
t in Fig. 3(a) for the same three values of a. The spe- 
cific heat jump at Tc, ACe(T'c)/7nTc is plotted versus a in 
Fig. 3(b). Figures 3(a) and 3(b) show that ACc{Tc)/-/nTc 
increases rather strongly with increasing a. 



Since A(i) and therefore dA^/(ii|t=i are assumed to be 
the same in the a-model as in the BCS theory, Eqs. (21) 
and (22) give the heat capacity jump at Tc as 



ACc(Tc) ^ ACc(Tc) 

JnTc 



12 
1.426 



BCS 

a 

aecs 

a 

"BCS 



CkBCS 



(50) 



The proportionality ACc(Tc)/7nTc ex a^ was previously 
noted. ^'* The numerical calculations of Ccs{t)/t^nTc at 
t = \ for the a-model for various values of a in Fig. 1 of 
Ref. 14 also agree with Eq. (50). 

The T = value ifc(O) of the thermodynamic critical 
field is given by Eq. (30) as 



gc(0 ) 



1/2 



1.382 a. 



(51) 



The dependences of Hc/Hc{0) on t and t^ obtained by 
replacing aecs by a in numerical solutions of Eqs. (27) 
and (29) are plotted in Figs. 4(a) and 4(b), respectively, 
for five a values as listed including aecs • On increasing 
a through a w 2. the curvature of Hc{t) / Hc{0) versus t"^ 
in Fig. 4(b) changes from positive to negative and the de- 
viation from a 1 — f^ dependence from negative to positive 
for < f < 1, as previously documented.^ 

Changing a from the BCS value to a different value 
involves more than just multiplying the BCS thermody- 
namic results by a power of a/aecs j since in addition to 
its presence in a prefactor, the parameter a is embedded 
within the respective integrals via the Fermi function in 
Eq. (20c). Therefore, e.g., one cannot use a table of BCS 
thermodynamic property values versus T to calculate the 
corresponding T-dependent predictions of the a-model. 
For example: (1) This is clear from Figs. 2, 3(a) and 4 
where the shapes of the respective curves versus tempera- 
ture strongly depend on a; (2) The superconducting and 
normal state entropies at Tc are not the same if one sim- 
ply multiplies the BCS superconducting state S'se(i)/7nTc 
values in Eq. (20a) or Fig. 2 by (a/ascs)^ ¥" li which is 
in conflict with our numerical calculations above which 
demonstrate that the a-model predicts a second-order 
transition at Tc; and (3) From Eq. (20b), if one ignores 
the presence of a inside the integral, one would infer that 
ACc(Tc)/7nTc ex a^ instead of the exact prediction with 
an a^ dependence in Eq. (50). 

Using Eq. (50), an accurate estimate of a can be ob- 
tained within the a-model from an accurate measurement 
of a sharp bulk heat capacity jump at Tc. However, the 
above discussion shows that to obtain accurate predic- 
tions of the entropy, heat capacity and thermodynamic 
critical field versus temperature for a given value of a, 
one should do numerical calculations for that a using the 
appropriate equations. Values of S'cs/7nTc, Ccs/jnTc and 
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Hc/Hc{0) versus t are given in the Appendix for seven a 
values (including the BCS predictions with a = aecs) in 
Tables III, IV and V, respectively. 

Perhaps surprisingly, if aecs is consistently changed 
to be a fixed but arbitrary value a throughout the calcu- 
lations, including those for A{t) and dA'^{t)/dt, we find 
that Scs{t)/jnTc and Ccs(i)/7nT'c are independent of a. 
In other words, if the a- model is treated self-consistently, 
the thermodynamic properties are the same as already 
presented for the BCS theory in Sees. IV and V for 
ct = aBCS ~ 1.764, irrespective of the value of a. 



VIII. ORIGINS OF DEVIATIONS OF THE 

HEAT CAPACITY JUMP AT Tc FROM THE BCS 

PREDICTION 



TABLE I: Heat capacity jump at Tc, ACc(rc), normalized 
by the BCS value ACc(rc)BCS for an isotropic s-wave gap 
in Eq. (22), and for superconductors with anisotropic gaps 
A = AoF{9, (f)) as in Eq. (52), where Ao is a positive constant 
and F{6, (p) is its angular dependence normalized such that 
the second moment is {F'^{d,(f>))FS = 1- The dimensionality 
of the FS is 2D (cyhndrical) and/or 3D (spherical). 



order parameter ^"(^1 



ACe(rc)/ACe(Tc)BCS 



2D, 3D s-wave 
2D anis. s-wave 

2D d-wave 

3D axial s-wave 



1 

a + cos^(2(j>) 
-^(3+8a + 8a2)/8 

(a<l) 
V2 cos(20) 



1 

2[3 + 8a-|-8n^]^ 



yiT^T? 



35-|-32a(l + a){5-|-4a+4a2) 

2/3 

7(15 + 10a + 3a^)^ 

5(315+420a+378a^ + 180a^+35a*) 



In this section we discuss superconducting gap 
anisotropics and other effects that can give rise to dif- 
ferences in the thermodynamic properties of supercon- 
ductors from the predictions of the BCS theory that 
the single-band a-model is generically formulated to fit. 
Thus from the fits one can quantify such deviations, 
which can then be interpreted in terms of other models. 

If the superconducting gap is anisotropic in mo- 
mentum space or scattering of the conduction carri- 
ers by magnetic impurities occurs, then one can ob- 
tain a < aBCS and ACes{Tc)/^nTc < 1.43 in a single- 
band superconductor. ^'^^ Openov has presented a com- 
prehensive theory on the effect of anisotropy in the su- 
perconducting order parameter on the specific heat jump 
ACc/jnTc of BCS superconductors (i.e., within the weak- 
coupling mean-field approximation).^ The calculations 
were carried out not only for clean superconductors as in 
the original BCS theory but also for those in which non- 
magnetic and/or magnetic impurity scattering of the con- 
duction electrons occurs. We mainly discuss the clean- 
limit predictions here. 

The wave vector k dependent order parameter Ak is 
defined as 



Ak = AoF{9, 



(52) 



where Ao(T) is an angle- independent positive constant 
and F gives its angular dependence on the Fermi surface 
(FS) in cyhndrical [2D, F = F{(j))] or spherical [3D, F = 
F{9, (p)] coordinates where 9 and (/> arc the polar and 
azimuthal angles, respectively. Then the prediction is^ 



ACe(rc) 

JnTc 



12 



7C(3) 



:f'{9,^))1, 
[FHd,' 



;fs 



(53) 



where (• • • )fs is an average of the enclosed quantity over 
the Fermi surface and the prefactor in square brackets 
is the heat capacity jump ACc{Tc)Bcs/lnTc ~ 1.426 in 
Eq. (22) for a BCS superconductor with an isotropic gap 



(F = 1). The moments of F are defined as 

(i^"(0))Fs = ^y'"d0i^"(0) (2D), (54) 

(F"(0,0))fs = ^j\9sxn9J\<pF%e,cjy) (3D). 

We normalize F{9,(t)) such that {F'^{6,(f)))Fs = 1, and 
then rewrite Eq. (53) as 



ACc(Tc) 

ACe(Te)BCS { FH 



;fs 



(55) 



For simplicity, we consider here free-electron Fermi sur- 
faces in either two dimensions (2D, cylindrical) or three 
(3D, spherical). Several types of gap anisotropy for 2D 
square or 3D tetragonal symmetry^^ arc considered as 
shown in Table I. Polar plots of the angular dependences 
of the 2D anisotropic s-wave and d-wave order parame- 
ters are shown in Fig. 7. Conventional superconductors 
driven by the electron-phonon interaction have s-wave 
order parameters with the same sign everywhere on the 
Fermi surface, whereas some unconventional supercon- 
ductors such as the high-Tc cuprates have sign-changing 
d-wave order parameters. ^^ An example of a 3D axial 
anisotropic s-wave gap^^ is shown in Fig. 8. 

The heat capacity jump at T^ in these cases is 
calculated using Eq. (55). For 2D d-wave pairing, 
ACc(rc)/ACc(rc)Bcs has the single value of 2/3 (Ta- 
ble I). Plots of ACc(rc)/ACc(Tc)Bcs versus the 
respectively-defined anisotropy parameter a for 2D 
anisotropic s-wave and 3D axial anisotropic s-wave order 
parameters calculated from the expressions in Table I us- 
ing Eq. (55) are shown in Fig. 9. One sees that the heat 
capacity jump at Tc monotonically decreases from the 
BCS value with increasing anisotropy, which corresponds 
to decreasing a for 2D anisotropic s-wave and increas- 
ing a for 3D axial anisotropic s-wave order parameters. 
The maximum decrease in the limit of large anisotropy 
is about a factor of two for both types of anisotropy. 
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a = 0.5 



anisotropic 




FIG. 7: (Color online) (top panel) Angular dependences of 
the 2D anisotropic s-wave order parameters with anisotropy 
parameters a = 0.5 (red) and a = 0.5 (blue) versus azimuthal 
angle cj> measured from the positive fc^-axis. (bottom panel) 
2D d-wave superconducting order parameter versus (j). The 
anisotropic s-wave order parameters have the same sign on all 
parts of the Fermi surface, whereas the d-wave order param- 
eter changes sign versus 4> ^^ shown. The data are computed 
from the expressions in Table I. 




FIG. 8: (Color online) Angular dependence of the 3D cylin- 
drically symmetric but axially anisotropic s-wave order pa- 
rameter in Table I for anisotropy parameter a = 2. The polar 
z-axis is vertical. 
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As mentioned in the introduction, values of the spe- 
cific heat jump ACo(T'c)/7n7c larger than the BCS value 
of 1.43 are observed in moderate- and strong-coupling 
superconductors for which a > aecs-^ Combescot con- 
cluded that a has the upper limit a < 2 within the weak- 
coupling BCS theory. ^'^ When magnetic impurities are 
present in a BCS superconductor, both Tc and the nor- 
malized specific heat jump ACc(Tc)/"fnTc are reduced 
compared to the BCS values with increasing electron- 
impurity scattering rate.^'^^ Both quantities go to zero 
at sufficiently high scattering rates. Nonmagnetic im- 
purites can also depress both Tc and ACc(Tc)/7nTc for 
single-band materials with anisotropic superconducting 
order parameters.^ 



FIG. 9: (Color online) Heat capacity jump ACe(Tc) at Tc nor- 
malized by the BCS value ACo(rc)BCS versus the anisotropy 
parameter a in Table I for the 2D anisotropic s-wave and the 
3D axial anisotropic s-wave order parameters. The different 
behaviors versus a result from the different roles of a in the 
anisotropy expressions in Table I. 



IX. EXTENSIONS OF THE SINGLE-BAND 
a-MODEL 



The single-band a-model in the clean limit has been 
extended to the two-band a-model in which each of two 
electron bands develop distinct isotropic superconducting 
energy gaps on the two respective Fermi surfaces. The 
parameter a = A{0) / kBTc and A(0) are different for the 
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two bands but Tc is the same. Within the weak-coupUng 
theory, Eq. (53) for the heat capacity jump at T^ for an 
anisotropic single gap is generalized to^° 



model to fit the temperature dependences of the Lon- 
don penetration depths of muhi-gap FeAs-based high-Tc 
superconductors.^^ 



AC^c(rc) 


12 

i7C(3)J 


[c,{Fl{0,^))ps+C2{Fiie,(j)))FS? 


7nTc 


c,{Fl{e,(b))FS+C2{Fi{e,cj)))Fs 

(56) 



where ci and C2 are the respective fractions of the total 
A^(0) contributed by bands 1 and 2, with ci + C2 = 1, 
and tlie prefactor in square brackets is the BCS value of 
ACe(rc)/7„Tc in Eq. (22). 

The isotropic s-wave two-band (two-gap) a-model has 
been applied to understand the temperature dependences 
of the specific heat and London penetration depth of ma- 
terials such as the fiducial two-gap compound MgB2 with 
Tc = 39 K (Refs. 14, 21, 22) and NbSe2 (Tc = 7.1 K).23 
It seems likely that many previously-studied supercon- 
ductors in the clean limit have multi-gap behaviors in 
addition to those already identified. Dolgov et al. made 
a detailed comparison of the predictions of the two- 
band a-modcl with the corresponding solution to the full 
Eliashberg equations and found good agreement. ^^ When 
strong interband scattering by nonmagnetic impurities is 
present, within the Eliashberg formalism the two gaps 
merge into a single gap.^^ Interband scattering by non- 
magnetic impurities has a similar strong effect on the 
superconductivity in a two-band, two-gap superconduc- 
tor as does intraband scattering by magnetic impurities 
in a single-gap superconductor. 

As noted above, the a-model is not self-consistent 
because a = aecs is used to calculate A{t) and the 
London penetration depth, but a different a value is 
used to calculate the thermodynamic properties. Ko- 
gan, Martin and Prozorov formulated a two-band (two- 
gap) weak-coupling clean-limit model with isotropic s- 
wave gaps on each Fermi surface, termed the 7-model, 
that self-consistently treats the temperature dependences 
of the two gaps and the associated London penetration 
depth and thermodymic properties. ^^ They find that 
ACo(Tc)/fcBT'c has an upper limit given by the BCS 
value of 1.43, and can be strongly suppressed from that 
value depending on the relative values of A^(0) of the 
two bands and the intraband and interband electron- 
phonon coupling constants. The same effect was doc- 
umented in Sec. VIII above for anisotropy of the gap 
within the single-band model. Thus when considering 
the effects of gap anisotropy on the heat capacity jump 
at Tc, one should evidently consider gap anisotropy glob- 
ally within the Brillouin zone. That is, when two Fermi 
surfaces each have isotropic gaps, but where the gaps 
are different on each Fermi surface, then the gap is 
wave-vector-dependent within the Brillouin zone, lead- 
ing to a suppression of the heat capacity jump from the 
BCS value qualitatively similar to that in the anisotropic 
single-band model. Kogan, Martin and Prozorov applied 
their 7-model to precisely fit the temperature depen- 
dences of the London penetration depth and electronic 
heat capacity of MgB2 and VsSi (Tc — 17 K). Prozorov 
and Kogan have also reviewed applications of the same 
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SUMMARY 



The single-band a-model has been extensively used by 
experimentalists in the past to fit electronic heat capacity 
versus temperature data for superconductors that deviate 
from the BCS prediction. The model is based on the BCS 
theory and assumes the same values of the normalized 
gap and London penetration depth as predicted by BCS 
that are obtained using the BCS value q;bcs ~ 1.764 of 
a = A(0)/A:bTc. However, in calculating the electronic 
entropy, heat capacity, free energy and thermodynamic 
critical field of the superconducting state, the a-model 
takes a to be a variable, which then allows calculations 
of these thermodynamic quantities to be adjusted to fit 
experimental data. This is an inconsistency in the model. 

Most previous experimental papers fitting the a model 
to experimental thermodynamics data do not explain 
how the theoretical values used in the fit were obtained. 
We have written the equations for the thermodynamic 
predictions of the BCS theory in terms of a to clarify 
how to calculate these quantities versus a and T and 
have compared the results for different a with the BCS 
predictions. Tables of values of the BCS predictions of 
the superconducting order parameter, London parame- 
ter and London penetration depth versus temperature 
are given in the Appendix, which are the same in the 
a-model by assumption, and of the Pippard penetration 
depth versus temperature. Additional tables in the Ap- 
pendix give the temperature dependences of the super- 
conducting state electronic entropy, heat capacity and 
thermodynamic critical field for seven representative val- 
ues of a including aecs to facilitate fitting of experimen- 
tal data by the a-model. These results could be interpo- 
lated to obtain values for other values of a. These tables 
supplement the table in Ref. 8 of superconducting state 
properties versus temperature predicted by the BCS the- 
ory. 

We find that if the a-model is treated self-consistently, 
i.e., if the same value of a is used to calculate the temper- 
ature dependence of the superconducting gap and of the 
thermodynamic properties, then the normalized thermo- 
dynamic properties versus temperature are independent 
of a and are therefore the same as presented for the BCS 
theory in Sees. IV and V and Tables III-V in the Ap- 
pendix for a = aecs ~ 1.764. 

Mechanisms for producing deviations of the super- 
conducting state thermodynamic properties from the 
BCS predictions were discussed that can be quanti- 
fied using the a-model. It is well-known that strong 
electron-phonon coupling increases the heat capacity 
jump ACo(Tc)/7nTc at Tc compared to the BCS value 
of 1.43.^ We calculated the influence of superconducting 
gap anisotropy in momentum space on the specific heat 
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jump at Tc for several types of gap anisotropy using the 
formalism discussed in Ref. 6, and found that the jump 
monotonically decreases with increasing anisotropy from 
the BCS value by up to a factor of two. Extensions of 
the a model were also discussed, including the two-band 
a-model^" and the self-consistent two-band 7-model.^^'^^ 
The present work points the way towards a uniform 
application of the a-model by experimentalists in their 
analyses of superconducting state thermodynamic data. 
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Appendix: Tables of Values 

Table II gives the T dependence of the BCS supercon- 
ducting gap with high resolution in t, especially near Tc, 
a dependence that is retained in the a-model. The sub- 
sequent three tables give the T dependences of the elec- 
tronic entropy (Table III), electronic heat capacity (Ta- 
ble IV) and thermodynamic critical field (Table V) for 
seven values of a including the BCS value aecs ~ 1.764, 
where the exact value of aecs in Eq. (4) was used in com- 
puting the values in the tables. Values of the BCS Lon- 
don parameter and the London penetration depth versus 
t are given in Table VI and values of the Pippard pene- 
tration depth versus t are given in Table VII. Table VI 
also applies to the a-model. 
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TABLE II: Normalized superconducting order parameter A(t) 
within the BCS weak-couphng theory using Eq. (11). The data 


= A(f)/A(0) versus 
are the same in the 


reduced temperature t - 
a-model. 


= T/Tc calculated 


t 


A(f) 


t 


A(t) 


t 


A(t) 





1 


0.760000 


0.76399 


0.997246 


0.091038 


0.020000 


1.00000 


0.780000 


0.73863 


0.997545 


0.085956 


0.022763 


1.00000 


0.781224 


0.73701 


0.997812 


0.081157 


0.040000 


1.00000 


0.800000 


0.71104 


0.998050 


0.076625 


0.060000 


1.00000 


0.805016 


0.70374 


0.998262 


0.072345 


0.080000 


1.00000 


0.820000 


0.68095 


0.998451 


0.068303 


0.100000 


1.00000 


0.826220 


0.67103 


0.998620 


0.064487 


0.120000 


1.00000 


0.840000 


0.64801 


0.998770 


0.060883 


0.129036 


1.00000 


0.845118 


0.63906 


0.998904 


0.057480 


0.140000 


1.00000 


0.860000 


0.61173 


0.999023 


0.054268 


0.160000 


0.99999 


0.861962 


0.60797 


0.999129 


0.051234 


0.180000 


0.99996 


0.876973 


0.57786 


0.999224 


0.048370 


0.200000 


0.99988 


0.880000 


0.57148 


0.999308 


0.045666 


0.220000 


0.99971 


0.890352 


0.54880 


0.999383 


0.043113 


0.223753 


0.99967 


0.900000 


0.52634 


0.999450 


0.040702 


0.240000 


0.99941 


0.902276 


0.52084 


0.999510 


0.038426 


0.260000 


0.99892 


0.912904 


0.49400 


0.999563 


0.036277 


0.280000 


0.99818 


0.920000 


0.47491 


0.999611 


0.034249 


0.300000 


0.99712 


0.922375 


0.46829 


0.999653 


0.032334 


0.308169 


0.99659 


0.930817 


0.44371 


0.999691 


0.030525 


0.320000 


0.99569 


0.938340 


0.42024 


0.999725 


0.028818 


0.340000 


0.99382 


0.940000 


0.41484 


0.999755 


0.027206 


0.360000 


0.99146 


0.945046 


0.39787 


0.999781 


0.025685 


0.380000 


0.98854 


0.951022 


0.37657 


0.999805 


0.024248 


0.383405 


0.98799 


0.956348 


0.35631 


0.999826 


0.022892 


0.400000 


0.98503 


0.960000 


0.34161 


0.999845 


0.021612 


0.420000 


0.98088 


0.961095 


0.33705 


0.999862 


0.020403 


0.440000 


0.97603 


0.965326 


0.31877 


0.999877 


0.019262 


0.450459 


0.97320 


0.969097 


0.30141 


0.999890 


0.018184 


0.460000 


0.97044 


0.972458 


0.28495 


0.999902 


0.017167 


0.480000 


0.96407 


0.975453 


0.26935 


0.999913 


0.016207 


0.500000 


0.95688 


0.978122 


0.25456 


0.999922 


0.015300 


0.510221 


0.95288 


0.980000 


0.24358 


0.999931 


0.014444 


0.520000 


0.94883 


0.980502 


0.24056 


0.999938 


0.013637 


0.540000 


0.93986 


0.982622 


0.22730 


0.999945 


0.012874 


0.560000 


0.92993 


0.984512 


0.21476 


0.999951 


0.012154 


0.563484 


0.92810 


0.986196 


0.20288 


0.999956 


0.011474 


0.580000 


0.91899 


0.987697 


0.19165 


0.999961 


0.010832 


0.600000 


0.90699 


0.989035 


0.18103 


0.999965 


0.010226 


0.610955 


0.89995 


0.990228 


0.17099 


0.999969 


0.0096540 


0.620000 


0.89387 


0.991290 


0.16150 


0.999972 


0.0091140 


0.640000 


0.87957 


0.992238 


0.15252 


0.999975 


0.0086042 


0.653263 


0.86939 


0.993082 


0.14404 


0.999978 


0.0081229 


0.660000 


0.86401 


0.993834 


0.13602 


0.999981 


0.0076685 


0.680000 


0.84710 


0.994505 


0.12845 


0.999983 


0.0072395 


0.690970 


0.83723 


0.995102 


0.12129 


0.999985 


0.0068346 


0.700000 


0.82877 


0.995635 


0.11453 


0.999986 


0.0064523 


0.720000 


0.80890 


0.996110 


0.10815 


0.999988 


0.0060913 


0.724577 


0.80412 


0.996533 


0.10212 


0.999989 


0.0057506 


0.740000 


0.78735 


0.996910 


0.096419 


1 





0.754529 


0.77057 
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TABLE III: Normalized superconducting state electronic 
the Qf-model in Eq. (20a) for the q values listed, where a 



entropy Scs/ThTc versus 
= A(0)/fcBT. The exact 



reduced temperature t — T/Tc predicted by 
value of qbcs ~ 1.764 is given in Eq. (4). 











Scb/IuTc 








t 


a = 1 


a = 1.25 


Q = 1.5 


a = aacs 


a = 2 


a = 2.25 


a = 2.5 


0.00 























0.02 


1.0784e-21 


5.5751e-27 


2.7177e-32 


6.4290e-38 


5.77506-43 


2.5627e-48 


1.11676-53 


0.04 


5.6943e-ll 


1.5142e-13 


3.8053e-16 


6.5731e-19 


2.15686-21 


4.9478e-24 


1.11506-26 


0.06 


2.0045e-07 


4.2516e-09 


8.5414e-ll 


1.3255e-12 


3.10426-14 


5.7079e-16 


1.03146-17 


0.08 


1.1596e-05 


6.9241e-07 


3.9244e-08 


1.8224e-09 


1.13896-10 


5.9230e-12 


3.02846-13 


0.10 


1.3076e-04 


1.4495e-05 


1.5281e-06 


1.3680e-07 


1.53916-08 


1.4926e-09 


1.42366-10 


0.12 


6.5353e-04 


1.0923e-04 


1.7394e-05 


2.4094e-06 


4.00916-07 


5.8866e-08 


8.50336-09 


0.14 


2.0559e-03 


4.6012e-04 


9.8269e-05 


1.8577e-05 


4.08566-06 


8.0636e-07 


1.56626-07 


0.16 


4.8480e-03 


1.3493e-03 


3.5885e-04 


8.5599e-05 


2.31956-05 


5.7125e-06 


1.38516-06 


0.18 


9.4393e-03 


3.1107e-03 


9.8060e-04 


2.8013e-04 


8.92516-05 


2.6103e-05 


7.51856-06 


0.20 


1.6080e-02 


6.0633e-03 


2.1889e-03 


7.2207e-04 


2.61806-04 


8.7833e-05 


2.90296-05 


0.22 


2.4863e-02 


1.0465e-02 


4.2201e-03 


1.5656e-03 


6.30836-04 


2.3676e-04 


8.75606-05 


0.24 


3.5758e-02 


1.6494e-02 


7.2929e-03 


2.9834e-03 


1.31256-03 


5.4078e-04 


2.19636-04 


0.26 


4.8656e-02 


2.4247e-02 


1.1590e-02 


5.1499e-03 


2.44036-03 


1.0881e-03 


4.78356-04 


0.28 


6.3395e-02 


3.3758e-02 


1.7251e-02 


8.2286e-03 


4.15596-03 


1.9830e-03 


9.33076-04 


0.30 


7.9798e-02 


4.5005e-02 


2.4372e-02 


1.23646-02 


6.60006-03 


3.3400e-03 


1.66726-03 


0.32 


9.7683e-02 


5.7933e-02 


3.3012e-02 


1.7678e-02 


9.90646-03 


5.2791e-03 


2.77536-03 


0.34 


1.1688e-01 


7.2463e-02 


4.3197e-02 


2.4267e-02 


1.41986-02 


7.9208e-03 


4.36006-03 


0.36 


1.3722e-01 


8.8503e-02 


5.4927e-02 


3.2209e-02 


1.95856-02 


1.1383e-02 


6.52866-03 


0.38 


1.5857e-01 


1.0596e-01 


6.8184e-02 


4.1557e-02 


2.61636-02 


1.5778e-02 


9.39146-03 


0.40 


1.8079e-01 


1.2472e-01 


8.2937e-02 


5.2352e-02 


3.40166-02 


2.1214e-02 


1.30596-02 


0.42 


2.0379e-01 


1.4471e-01 


9.9144e-02 


6.4618e-02 


4.32146-02 


2.7790e-02 


1.76426-02 


0.44 


2.2746e-01 


1.6582e-01 


1.1676e-01 


7.83676-02 


5.38186-02 


3.5599e-02 


2.32496-02 


0.46 


2.5171e-01 


1.8798e-01 


1.3573e-01 


9.36036-02 


6.58766-02 


4.4728e-02 


2.99866-02 


0.48 


2.7649e-01 


2.1110e-01 


1.5600e-01 


1.10326-01 


7.94316-02 


5.5256e-02 


3.79576-02 


0.50 


3.0173e-01 


2.3510e-01 


1.7752e-01 


1.28526-01 


9.45176-02 


6.7257e-02 


4.72666-02 


0.52 


3.2736e-01 


2.5993e-01 


2.0024e-01 


1.48176-01 


1.11166-01 


8.0801e-02 


5.80116-02 


0.54 


3.5336e-01 


2.8552e-01 


2.2411e-01 


1.69266-01 


1.29386-01 


9.5951e-02 


7.02916-02 


0.56 


3.7968e-01 


3.1181e-01 


2.4907e-01 


1.91786-01 


1.49216-01 


1.1277e-01 


8.42026-02 


0.58 


4.0628e-01 


3.3876e-01 


2.7508e-01 


2.15696-01 


1.70656-01 


1.3131e-01 


9.98406-02 


0.60 


4.3314e-01 


3.6631e-01 


3.0209e-01 


2.40996-01 


1.93716-01 


1.5162e-01 


1.17306-01 


0.62 


4.6023e-01 


3.9443e-01 


3.3006e-01 


2.67636-01 


2.18406-01 


1.7377e-01 


1.36676-01 


0.64 


4.8753e-01 


4.2307e-01 


3.5895e-01 


2.95616-01 


2.44746-01 


1.9779e-01 


1.58056-01 


0.66 


5.1501e-01 


4.5221e-01 


3.8871e-01 


3.24886-01 


2.72716-01 


2.2373e-01 


1.81536-01 


0.68 


5.4267e-01 


4.8180e-01 


4.1932e-01 


3.55446-01 


3.02346-01 


2.5164e-01 


2.07206-01 


0.70 


5.7048e-01 


5.1183e-01 


4.5073e-01 


3.87266-01 


3.33616-01 


2.8156e-01 


2.35166-01 


0.72 


5.9843e-01 


5.4226e-01 


4.8292e-01 


4.20306-01 


3.66546-01 


3.1354e-01 


2.65516-01 


0.74 


6.2652e-01 


5.7306e-01 


5.1585e-01 


4.54566-01 


4.01116-01 


3.4761e-01 


2.98346-01 


0.76 


6.5473e-01 


6.0423e-01 


5.4949e-01 


4.89996-01 


4.37346-01 


3.8382e-01 


3.33756-01 


0.78 


6.8304e-01 


6.3573e-01 


5.8381e-01 


5.26596-01 


4.75216-01 


4.2221e-01 


3.71856-01 


0.80 


7.1146e-01 


6.6755e-01 


6.1880e-01 


5.64346-01 


5.14736-01 


4.6282e-01 


4.12736-01 


0.82 


7.3998e-01 


6.9967e-01 


6.5442e-01 


6.03206-01 


5.55896-01 


5.0568e-01 


4.56506-01 


0.84 


7.6858e-01 


7.3208e-01 


6.9064e-01 


6.43156-01 


5.98706-01 


5.5085e-01 


5.03276-01 


0.86 


7.9727e-01 


7.6475e-01 


7.2746e-01 


6.84196-01 


6.43156-01 


5.9836e-01 


5.53166-01 


0.88 


8.2603e-01 


7.9769e-01 


7.6484e-01 


7.26286-01 


6.89236-01 


6.4826e-01 


6.06286-01 


0.90 


8.5487e-01 


8.3086e-01 


8.0277e-01 


7.69416-01 


7.36956-01 


7.0057e-01 


6.62766-01 


0.92 


8.8377e-01 


8.6427e-01 


8.4123e-01 


8.13566-01 


7.86316-01 


7.5536e-01 


7.22716-01 


0.94 


9.1274e-01 


8.9789e-01 


8.8020e-01 


8.58716-01 


8.37296-01 


8.1265e-01 


7.86276-01 


0.96 


9.4177e-01 


9.3173e-01 


9.1966e-01 


9.04846-01 


8.89906-01 


8.7249e-01 


8.53586-01 


0.98 


9.7086e-01 


9.6577e-01 


9.5960e-01 


9.51956-01 


9.44146-01 


9.3492e-01 


9.24776-01 


1.00 


1 


1 


1 


1 


1 


1 


1 
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TABLE IV: Normalized sup 


erconducting 


state electronic heat 


capacity Ccb/JuTc 


versus reduced temperature t - 


= T/Tc predicted 


by the a 


-model in Eq. (20b; 


) for the a values listed, where a - 


= A(0)/fcBT. The exact value of 


ascs ~ 1.764 is 


given in Eq. (4). 










Cos/7nrc 








t 


a = 1 


a = 1.25 


a = 1.5 


a = ctBCS 


a = 2 


a = 2.25 


a = 2.5 


0.00 























0.02 


5.3420e-20 


3.4582e-25 


2.0254e-30 


5.6392e-36 


5.7472e-41 


2.8707e-46 


1.3904e-51 


0.04 


1.3992e-09 


4.6649e-12 


1.4098e-14 


2.8684e-17 


1.0684e-19 


2.7600e-22 


6.9165e-25 


0.06 


3.2619e-06 


8.6815e-08 


2.0988e-09 


3.8387e-ll 


1.0209e-12 


2.1147e-14 


4.2505e-16 


0.08 


1.4076e-04 


1.05518-05 


7.1992e-07 


3.9418e-08 


2.7984e-09 


1.6400e-10 


9.3298e-12 


0.10 


1.2644e-03 


1.7595e-04 


2.2336e-05 


2.3583e-06 


3.0150e-07 


3.2957e-08 


3.4980e-09 


0.12 


5.2492e-03 


l.lOlle-03 


2.1114e-04 


3.4498e-05 


6.5239e-06 


1.0799e-06 


1.7363e-07 


0.14 


1.4121e-02 


3.9645e-03 


1.0194e-03 


2.2733e-04 


5.6825e-05 


1.2645e-05 


2.7342e-06 


0.16 


2.9088e-02 


1.0151e-02 


3.2496e-03 


9.1427e-04 


2.8158e-04 


7.8195e-05 


2.1108e-05 


0.18 


5.0293e-02 


2.0771e-02 


7.8790e-03 


2.6544e-03 


9.6117e-04 


3.1697e-04 


1.0165e-04 


0.20 


7.7077e-02 


3.6407e-02 


1.5811e-02 


6.1494e-03 


2.5338e-03 


9.5848e-04 


3.5270e-04 


0.22 


1.0838e-01 


5.7117e-02 


2.7699e-02 


1.2113e-02 


5.5462e-03 


2.3469e-03 


9.6633e-04 


0.24 


1.4304e-01 


8.2570e-02 


4.3895e-02 


2.1163e-02 


1.0578e-02 


4.9137e-03 


2.2218e-03 


0.26 


1.8001e-01 


1.1221e-01 


6.4473e-02 


3.3758e-02 


1.8173e-02 


9.1348e-03 


4.4707e-03 


0.28 


2.1843e-01 


1.4540e-01 


8.9299e-02 


5.0188e-02 


2.8795e-02 


1.5487e-02 


8.1128e-03 


0.30 


2.5761e-01 


1.8149e-01 


1.1811e-01 


7.0594e-02 


4.2806e-02 


2.4419e-02 


1.3569e-02 


0.32 


2.9707e-01 


2.1990e-01 


1.5056e-01 


9.4995e-02 


6.0474e-02 


3.6326e-02 


2.1261e-02 


0.34 


3.3648e-01 


2.6012e-01 


1.8630e-01 


1.2333e-01 


8.1975e-02 


5.1554e-02 


3.1593e-02 


0.36 


3.7561e-01 


3.0174e-01 


2.2496e-01 


1.5547e-01 


1.0742e-01 


7.0389e-02 


4.4951e-02 


0.38 


4.1433e-01 


3.4441e-01 


2.6622e-01 


1.9127e-01 


1.3685e-01 


9.3070e-02 


6.1690e-02 


0.40 


4.5256e-01 


3.8785e-01 


3.0975e-01 


2.3054e-01 


1.7030e-01 


1.1980e-01 


8.2141e-02 


0.42 


4.9028e-01 


4.3186e-01 


3.5530e-01 


2.7313e-01 


2.0773e-01 


1.5073e-01 


1.0661e-01 


0.44 


5.2747e-01 


4.7626e-01 


4.0262e-01 


3.1884e-01 


2.4911e-01 


1.8601e-01 


1.3539e-01 


0.46 


5.6416e-01 


5.2093e-01 


4.5150e-01 


3.6752e-01 


2.9441e-01 


2.2575e-01 


1.6875e-01 


0.48 


6.0037e-01 


5.6576e-01 


5.0176e-01 


4.1899e-01 


3.4356e-01 


2.7006e-01 


2.0696e-01 


0.50 


6.3613e-01 


6.1069e-01 


5.5325e-01 


4.7313e-01 


3.9651e-01 


3.1903e-01 


2.5026e-01 


0.52 


6.7147e-01 


6.5566e-01 


6.0582e-01 


5.2978e-01 


4.5320e-01 


3.7275e-01 


2.9891e-01 


0.54 


7.0643e-01 


7.0062e-01 


6.5937e-01 


5.8882e-01 


5.1358e-01 


4.3129e-01 


3.5315e-01 


0.56 


7.4104e-01 


7.4556e-01 


7.1380e-01 


6.5014e-01 


5.7759e-01 


4.9473e-01 


4.1324e-01 


0.58 


7.7533e-01 


7.90448-01 


7.6900e-01 


7.1363e-01 


6.4518e-01 


5.6316e-01 


4.7943e-01 


0.60 


8.0934e-01 


8.3525e-01 


8.2492e-01 


7.7919e-01 


7.1630e-01 


6.3665e-01 


5.5198e-01 


0.62 


8.4309e-01 


8.7999e-01 


8.8149e-01 


8.4674e-01 


7.9091e-01 


7.1528e-01 


6.3115e-01 


0.64 


8.7660e-01 


9.2465e-01 


9.3865e-01 


9.1618e-01 


8.6897e-01 


7.9914e-01 


7.1724e-01 


0.66 


9.0991e-01 


9.6922e-01 


9.9634e-01 


9.8745e-01 


9.5045e-01 


8.8832e-01 


8.1051e-01 


0.68 


9.4303e-01 


1.0137 


1.0545 


1.0605 


1.0353 


9.8289e-01 


9.1129e-01 


0.70 


9.7598e-01 


1.0581 


1.1132 


1.1352 


1.1235 


1.0830 


1.0199 


0.72 


1.0088 


1.1024 


1.1722 


1.2115 


1.2150 


1.1886 


1.1366 


0.74 


1.0414 


1.1467 


1.2317 


1.2894 


1.3098 


1.3000 


1.2618 


0.76 


1.0740 


1.1908 


1.2915 


1.3689 


1.4078 


1.4171 


1.3959 


0.78 


1.1064 


1.2349 


1.3517 


1.4498 


1.5092 


1.5402 


1.5393 


0.80 


1.1387 


1.2789 


1.4122 


1.5321 


1.6137 


1.6693 


1.6923 


0.82 


1.1710 


1.3229 


1.4730 


1.6159 


1.7214 


1.8045 


1.8554 


0.84 


1.2031 


1.3668 


1.5341 


1.7010 


1.8324 


1.9460 


2.0291 


0.86 


1.2352 


1.4106 


1.5954 


1.7874 


1.9465 


2.0939 


2.2139 


0.88 


1.2673 


1.4544 


1.6570 


1.8750 


2.0638 


2.2483 


2.4102 


0.90 


1.2992 


1.4981 


1.7189 


1.9639 


2.1842 


2.4094 


2.6186 


0.92 


1.3311 


1.5418 


1.7809 


2.0541 


2.3078 


2.5774 


2.8398 


0.94 


1.3630 


1.5855 


1.8432 


2.1454 


2.4345 


2.7523 


3.0743 


0.96 


1.3948 


1.6291 


1.9057 


2.2378 


2.5644 


2.9343 


3.3229 


0.98 


1.4266 


1.6727 


1.9684 


2.3314 


2.6974 


3.1237 


3.5861 


1.00 


1.4584 


1.7162 


2.0313 


2.4261 


2.8335 


3.3205 


3.8649 
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TABLE V: Thermodynamic critical field normalized to its zero temperature value, Hdt)/ HdQ), versus reduced temperature 
t = T/Tc predicted by the a-model in Eq. (29) for the a values listed, where a = A(0)/fcBr. The exact value of OfBCS ~ 1.764 
is given in Eq. (4) and the expression for Hc{0) versus a is given in Eq. (51). 











mit)/m{o) 








t 


a = 1 


a = 1.25 


a = 1.5 


a = Qbcs 


Q = 2 


a = 2.25 


a = 2.5 


0.00 


1 


1 


1 


1 


1 


1 


1 


0.02 


0.99868 


0.99916 


0.99942 


0.99958 


0.99967 


0.99974 


0.99979 


0.04 


0.99472 


0.99663 


0.99766 


0.99831 


0.99868 


0.99896 


0.99916 


0.06 


0.98809 


0.99239 


0.99472 


0.99619 


0.99704 


0.99766 


0.99810 


0.08 


0.97872 


0.98643 


0.99060 


0.99321 


0.99472 


0.99583 


0.99663 


0.10 


0.96659 


0.97872 


0.98527 


0.98937 


0.99174 


0.99348 


0.99472 


0.12 


0.95172 


0.96924 


0.97872 


0.98466 


0.98809 


0.99060 


0.99239 


0.14 


0.93431 


0.95798 


0.97094 


0.97906 


0.98375 


0.98718 


0.98963 


0.16 


0.91468 


0.94503 


0.96192 


0.97256 


0.97871 


0.98321 


0.98642 


0.18 


0.89326 


0.93050 


0.95167 


0.96514 


0.97296 


0.97868 


0.98275 


0.20 


0.87047 


0.91456 


0.94025 


0.95680 


0.96646 


0.97355 


0.97860 


0.22 


0.84673 


0.89739 


0.92769 


0.94752 


0.95919 


0.96777 


0.97391 


0.24 


0.82238 


0.87917 


0.91406 


0.93730 


0.95110 


0.96132 


0.96863 


0.26 


0.79768 


0.86008 


0.89943 


0.92613 


0.94218 


0.95413 


0.96271 


0.28 


0.77286 


0.84025 


0.88387 


0.91402 


0.93237 


0.94615 


0.95609 


0.30 


0.74806 


0.81983 


0.86744 


0.90098 


0.92167 


0.93733 


0.94870 


0.32 


0.72340 


0.79891 


0.85020 


0.88702 


0.91004 


0.92764 


0.94050 


0.34 


0.69895 


0.YYY58 


0.83221 


0.87216 


0.89748 


0.91702 


0.93142 


0.36 


0.67475 


0.75591 


0.81352 


0.85640 


0.88397 


0.90547 


0.92142 


0.38 


0.65082 


0.73396 


0.79418 


0.83979 


0.86951 


0.89294 


0.91047 


0.40 


0.62718 


0.71177 


0.77422 


0.82232 


0.85410 


0.87941 


0.89853 


0.42 


0.60384 


0.68937 


0.75368 


0.80403 


0.83774 


0.86488 


0.88556 


0.44 


0.58077 


0.66680 


0.73261 


0.78494 


0.82044 


0.84933 


0.87154 


0.46 


0.55798 


0.64408 


0.71102 


0.76506 


0.80221 


0.83276 


0.85646 


0.48 


0.53545 


0.62121 


0.68895 


0.74443 


0.78305 


0.81515 


0.84028 


0.50 


0.51316 


0.59823 


0.66642 


0.72306 


0.76298 


0.79650 


0.82299 


0.52 


0.49110 


0.57514 


0.64346 


0.70098 


0.74200 


0.77682 


0.80459 


0.54 


0.46926 


0.55195 


0.62008 


0.67820 


0.72014 


0.75610 


0.78505 


0.56 


0.44761 


0.52867 


0.59631 


0.65474 


0.69740 


0.73434 


0.76437 


0.58 


0.42615 


0.50531 


0.57216 


0.63063 


0.67380 


0.71156 


0.74254 


0.60 


0.40486 


0.48186 


0.54766 


0.60588 


0.64934 


0.68774 


0.71954 


0.62 


0.38373 


0.45834 


0.52281 


0.58051 


0.62405 


0.66290 


0.69537 


0.64 


0.36274 


0.43475 


0.49763 


0.55454 


0.59793 


0.63703 


0.67002 


0.66 


0.34188 


0.41108 


0.47213 


0.52798 


0.57101 


0.61014 


0.64347 


0.68 


0.32115 


0.38736 


0.44633 


0.50085 


0.54328 


0.58224 


0.61573 


0.70 


0.30053 


0.36356 


0.42024 


0.47317 


0.51477 


0.55333 


0.58678 


0.72 


0.28002 


0.33971 


0.39387 


0.44494 


0.48548 


0.52340 


0.55661 


0.74 


0.25960 


0.31579 


0.36723 


0.41619 


0.45543 


0.49248 


0.52522 


0.76 


0.23927 


0.29182 


0.34033 


0.38693 


0.42463 


0.46055 


0.49259 


0.78 


0.21901 


0.26779 


0.31317 


0.35717 


0.39308 


0.42762 


0.45871 


0.80 


0.19884 


0.24371 


0.28578 


0.32692 


0.36081 


0.39369 


0.42358 


0.82 


0.17873 


0.21957 


0.25814 


0.29619 


0.32781 


0.35877 


0.38717 


0.84 


0.15868 


0.19537 


0.23028 


0.26500 


0.29411 


0.32286 


0.34949 


0.86 


0.13869 


0.17113 


0.20220 


0.23336 


0.25970 


0.28596 


0.31050 


0.88 


0.11875 


0.14683 


0.17391 


0.20127 


0.22461 


0.24807 


0.27021 


0.90 


0.09886 


0.12248 


0.14541 


0.16875 


0.18883 


0.20919 


0.22859 


0.92 


0.07902 


0.09808 


0.11670 


0.13581 


0.15238 


0.16932 


0.18563 


0.94 


0.05921 


0.07363 


0.08781 


0.10246 


0.11526 


0.12847 


0.14131 


0.96 


0.03944 


0.04914 


0.05872 


0.06870 


0.07749 


0.08664 


0.09561 


0.98 


0.01971 


0.02459 


0.02945 


0.03454 


0.03906 


0.04381 


0.04851 


1.00 
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TABLE VI: 


Normalized London paramet 


er I{t) = 1 - ^ from Eq. 


(39) and normalized London penetration depth ^hy = 


Al(«) 1 _ 


1 1 c — ™ t:i„ ^/l/l„^ , 


rersus reduced temperature 


t = T/Tc 


according to the BCS theory. 


' 1 ' ri ci rid + ci ara ^"rio 


>'l(0) 

same in the 


a-model. 


X lit UcLLcL die Liic 


t 


I{t) 


AAL(i)/AL(0) 


t 


lit) 


AAL(t)/AL(0) 


0.06 


2.3513e-12 


1.1757e-12 


0.54 


2.1335e-01 


1.2748e-01 


0.07 


1.4543e-10 


7.2714e-ll 


0.55 


2.2597e-01 


1.3664e-01 


0.08 


3.1803e-09 


1.5902e-09 


0.56 


2.3889e-01 


1.4624e-01 


0.09 


3.4812e-08 


1.7406e-08 


0.57 


2.5209e-01 


1.5631e-01 


0.10 


2.3490e-07 


1.1745e-07 


0.58 


2.6557e-01 


1.6687e-01 


0.11 


1.1155e-06 


5.5774e-07 


0.59 


2.7932e-01 


1.7795e-01 


0.12 


4.0716e-06 


2.0358e-06 


0.60 


2.9333e-01 


1.8957e-01 


0.13 


1.2142e-05 


6.0713e-06 


0.61 


3.0759e-01 


2.0176e-01 


0.14 


3.0902e-05 


1.5451e-05 


0.62 


3.2211e-01 


2.1456e-01 


0.15 


6.9287e-05 


3.4645e-05 


0.63 


3.3686e-01 


2.2800e-01 


0.16 


1.4018e-04 


7.0099e-05 


0.64 


3.5185e-01 


2.4211e-01 


0.17 


2.6064e-04 


1.3035e-04 


0.65 


3.6706e-01 


2.5695e-01 


0.18 


4.5173e-04 


2.2594e-04 


0.66 


3.8249e-01 


2.7256e-01 


0.19 


7.3799e-04 


3.6920e-04 


0.67 


3.9814e-01 


2.8899e-01 


0.20 


1.1467e-03 


5.7386e-04 


0.68 


4.1399e-01 


3.0631e-01 


0.21 


1.7071e-03 


8.5464e-04 


0.69 


4.3004e-01 


3.2458e-01 


0.22 


2.4491e-03 


1.2268e-03 


0.70 


4.4629e-01 


3.4387e-01 


0.23 


3.4029e-03 


1.7058e-03 


0.71 


4.6272e-01 


3.6427e-01 


0.24 


4.5978e-03 


2.3069e-03 


0.72 


4.7934e-01 


3.8587e-01 


0.25 


6.0616e-03 


3.0447e-03 


0.73 


4.9613e-01 


4.0878e-01 


0.26 


7.8204e-03 


3.9333e-03 


0.74 


5.1310e-01 


4.3311e-01 


0.27 


9.8978e-03 


4.9859e-03 


0.75 


5.3023e-01 


4.5901e-01 


0.28 


1.2315e-02 


6.2149e-03 


0.76 


5.4753e-01 


4.8663e-01 


0.29 


1.5090e-02 


7.6316e-03 


0.77 


5.6498e-01 


5.1616e-01 


0.30 


1.8240e-02 


9.2465e-03 


0.78 


5.8258e-01 


5.4779e-01 


0.31 


2.1776e-02 


1.1069e-02 


0.79 


6.0033e-01 


5.8179e-01 


0.32 


2.5710e-02 


1.3109e-02 


0.80 


6.1822e-01 


6.1843e-01 


0.33 


3.0051e-02 


1.5373e-02 


0.81 


6.3625e-01 


6.5805e-01 


0.34 


3.4803e-02 


1.7870e-02 


0.82 


6.5441e-01 


7.0106e-01 


0.35 


3.9972e-02 


2.0606e-02 


0.83 


6.7270e-01 


7.4795e-01 


0.36 


4.5559e-02 


2.3589e-02 


0.84 


6.9112e-01 


7.9931e-01 


0.37 


5.1565e-02 


2.6824e-02 


0.85 


7.0966e-01 


8.5587e-01 


0.38 


5.7988e-02 


3.0319e-02 


0.86 


7.2832e-01 


9.1853e-01 


0.39 


6.4827e-02 


3.4080e-02 


0.87 


7.4709e-01 


9.8846e-01 


0.40 


7.2078e-02 


3.8112e-02 


0.88 


7.6597e-01 


1.0671 


0.41 


7.9736e-02 


4.2423e-02 


0.89 


7.8496e-01 


1.1565 


0.42 


8.7797e-02 


4.7018e-02 


0.90 


8.0405e-01 


1.2591 


0.43 


9.6255e-02 


5.1906e-02 


0.91 


8.2324e-01 


1.3786 


0.44 


1.0510e-01 


5.7093e-02 


0.92 


8.4253e-01 


1.5200 


0.45 


1.1433e-01 


6.2587e-02 


0.93 


8.6192e-01 


1.6911 


0.46 


1.2394e-01 


6.8397e-02 


0.94 


8.8139e-01 


1.9037 


0.47 


1.3391e-01 


7.4532e-02 


0.95 


9.0096e-01 


2.1775 


0.48 


1.4425e-01 


8.1000e-02 


0.96 


9.2061e-01 


2.5490 


0.49 


1.5493e-01 


8.7813e-02 


0.97 


9.4034e-01 


3.0940 


0.50 


1.6596e-01 


9.4982e-02 


0.98 


9.6015e-01 


4.0093 


0.51 


1.7733e-01 


1.0252e-01 


0.99 


9.8004e-01 


6.0775 


0.52 


1.8902e-01 


1.1044e-01 


1.00 


1 




0.53 


2.0103e-01 


1.1875e-01 
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TABLE VII: Normalized Pippard penetration depth ^ ^,L' 
from Eq. (46a) versus reduced temperature i = T/T^ accord- 
ing to the BCS theory. 



t 



AAp(t)/Ap(0) 



AAp(i)/Ap(0) 



0.16 


1.4595e-05 


0.58 


7.1394e-02 


0.17 


2.8455e-05 


0.59 


7.6631e-02 


0.18 


5.1314e-05 


0.60 


8.2150e-02 


0.19 


8.6852e-05 


0.61 


8.7964e-02 


0.20 


1.3942e-04 


0.62 


9.4090e-02 


0.21 


2.1399e-04 


0.63 


1.0054e-01 


0.22 


3.1601e-04 


0.64 


1.0734e-01 


0.23 


4.5135e-04 


0.65 


1.1451e-01 


0.24 


6.2615e-04 


0.66 


1.2207e-01 


0.25 


8.4673e-04 


0.67 


1.3005e-01 


0.26 


1.1195e-03 


0.68 


1.3847e-01 


0.27 


1.4509e-03 


0.69 


1.4736e-01 


0.28 


1.8472e-03 


0.70 


1.5676e-01 


0.29 


2.3148e-03 


0.71 


1.6671e-01 


0.30 


2.8596e-Q3 


0.72 


1.7724e-01 


0.31 


3.4876e-Q3 


0.73 


1.8842e-01 


0.32 


4.2046e-03 


0.74 


2.0028e-01 


0.33 


5.Q161e-Q3 


0.75 


2.1290e-01 


0.34 


5.9275e-Q3 


0.76 


2.2635e-01 


0.35 


6.9441e-03 


0.77 


2.4069e-01 


0.36 


8.Q7Q9e-03 


0.78 


2.5603e-01 


0.37 


9.3129e-03 


0.79 


2.7247e-01 


0.38 


1.0675e-02 


0.80 


2.9014e-01 


0.39 


1.2162e-Q2 


0.81 


3.0918e-01 


0.40 


1.3780e-02 


0.82 


3.2977e-01 


0.41 


1.5531e-02 


0.83 


3.5211e-01 


0.42 


1.7423e-Q2 


0.84 


3.7645e-01 


0.43 


1.9459e-02 


0.85 


4.0309e-01 


0.44 


2.1645e-02 


0.86 


4.3240e-01 


0.45 


2.3985e-Q2 


0.87 


4.6486e-01 


0.46 


2.6486e-02 


0.88 


5.0106e-01 


0.47 


2.9153e-Q2 


0.89 


5.4177e-01 


0.48 


3.1991e-Q2 


0.90 


5.8800e-01 


0.49 


3.5007e-02 


0.91 


6.4115e-01 


0.50 


3.82Q8e-Q2 


0.92 


7.0315e-01 


0.51 


4.16Qle-02 


0.93 


7.7683e-01 


0.52 


4.5192e-02 


0.94 


8.6653e-01 


0.53 


4.8989e-02 


0.95 


9.7930e-01 


0.54 


5.30Q2e-Q2 


0.96 


1.1277 


0.55 


5.7238e-02 


0.97 


1.3370 


0.56 


6.1709e-Q2 


0.98 


1.6697 


0.57 


6.6423e-02 


0.99 


2.3567 



